Staggered flux vortices and the superconducting transition in the layered cuprates 
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We propose an effective model for the superconducting transition in the high-T c cuprates moti- 
vated by the SU(2) gauge theory approach. In addition to variations of the superconducting phase 
we allow for local admixture of staggered flux order. This leads to an unbinding transition of vortices 
with staggered flux core that are energetically preferable to conventional vortices. Based on param- 
eter estimates for the two-dimensional t-J-model we argue that the staggered flux vortices provide 
a way to understand a phase with a moderate density of mobile vortices over a large temperature 
range above T c that yet exhibits otherwise normal transport properties. This picture is consistent 
with the large Nernst signal observed in this region. 



The nature of the pseudogap phase of the (hole-)- 
underdoped high-T c cuprates is one of the central ques- 
tions of correlated electron physics. A number of sce- 
narios and descriptions are successful in capturing cer- 
tain aspects of the problem. Yet many theories face sub- 
stantial difficulties when it comes to combining the large 
number of experimentally established anomalies of the 
underdoped state. A new challenge in this context has 
been set out recently by the Nernst effect measurements 
in underdoped samples^]. In these experiments a ther- 
mal gradient is applied in the copper-oxide planes. In the 
presence of a small out-of-plane magnetic field, a voltage 
drop in the direction perpendicular to magnetic field and 
thermal gradient is observed. This voltage is interpreted 
as the phase slip signal arising from vortices moving from 
hot to cold. Thus the Nernst effect experiment reveals 
the existence of substantial superconducting (SC) short 
range correlations over a sizable region, starting signifi- 
cantly below the onset temperature T* of the pseudogap, 
but extending up to temperatures high above the low T c s 
of underdoped samples. 

At first sight the observation of vortices above T c fits 
well into a scenario^ where SC phase fluctuations de- 
stroy the long range coherence, and short range pairing 
correlations survive up to much higher temperatures. In 
view of the small superfiuid weight of underdoped sys- 
tems it is conceivable that the SC transition is driven 
by a vortex unbinding similar to the XY-transition, the 
critical temperature disappearing like the doping x for 
x — > 0. However a simple phase fluctuation scenario faces 
the following problem. The creation of a vortex comes 
at a price, as the SC order parameter goes to zero in 
the vortex core and condensation energy is lost. In dis- 
ordered films of conventional superconductors, where a 
Berezinskii-Kosterlitz-Thouless(BKT) transition can be 
observed, the mean-field critical temperature and the 
vortex unbinding temperature are similar, and due to 
the mean free path I entering the effective coherence 
length £ = \/£$t (with £o ~ vp/A) the vortex core en- 
ergy ~ £ 2 A 2 /ef becomes small near T C) as A decreases 
faster for T — ► T c than £ increases. Thus the loss of 
condensation energy is limited and the vortices can be 
cheap. In contrast with that the high-T c cuprates are in 
the clean limit and for underdoped samples the gap mag- 



nitude remains large up to temperatures far above the 
SC transition. Hence the vortex core energy for bringing 
the gap magnitude down to zero inside the vortex would 
normally be expected to be huge (sf in BCS theory and 
of order J in our case), and only exponentially few of 
these expensive vortices could be created above the small 
T c in underdoped samples. Then we would expect that 
the transport properties of the pseudogap state resemble 
those of a flux-flow (FF) phase. Transport resembling 
the normal state only occurs at high temperatures of or- 
der of the core energy when the vortices proliferate and 
overlap. However this picture is inconsistent with ex- 
periments. These show that above a limited fluctuation 
regime close to T c the in-plane transport looks rather 
normal and signs of FF conductivity do not extend far 
above T c - contrary to the Nernst signal. Apparently the 
conductivity a n due to a significant number of quasiparti- 
cles dominates over the FF conductivity <7ff in the total 
conductivity a = a n + <tff • Thus one has to explain two 
things: where the normal excitations come from and why 
the FF contribution is small, off can be estimated to be 
oc rj/ny, where ny is the density of vortices either forced 
in by a magnetic field in the mixed state or generated 
thermally above the BKT transition. n is the friction co- 
efficient for the vortex motion. The FF conductivity is 
small if 1] is small and ny is not. In the following we show 
that our model produces a moderately large ny even for 
underdoped samples with low T c . We will also present 
an approximate calculation that yields a finite density of 
normal excitations. We will not attempt to calculate n in 
this work. An effect that may reduce n is the experimen- 
tally observed ;3] small low-energy density of states in the 
vortex cores. This translates into small dissipation due 
to vortex motion and thus small off- Note that n can- 
not become arbitrarily small if we want a n to dominate 
the conductivity, as the the quasiparticles responsible for 
<j n will cause some dissipation. Thus a small FF contri- 
bution requires a moderately large ny at temperatures 
above the limited fluctuation regime near T c . 

An extreme way to obtain normal transport properties 
just above T c is to make the vortices very dense, such that 
they overlap just above the transition. Then it is hard 
to understand why the vortex Nernst signal persists to 
temperatures so high above T c . A theory considering 
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purely Gaussian SC fluctuations gave good agreement 
for the Nernst effect in overdoped and optimally doped 
samples. However the description of underdoped samples 
becomes problematic and an additional suppression of T c 
had to be invoked. 

Thus neither very few and expensive nor too many and 
too cheap vortices seem to match the experimental pic- 
ture. Basically what is needed is a theory which produces 
a core energy of order T c rather than J. In other words we 
need to decrease the core energy by placing in the vortex 
core another non-superconducting state that is nearby in 
energy. There are several proposals in the literature]^ 
for such a cheap vortex core, mainly emphasizing the 
vicinity of the d-wave superconductor to other ordered 
states. Here we study the possibility of a staggered flux 
(SF) state inside the vortexja, Qj ■ This scenario has the 
advantage that it emerges naturally from the SU(2) in- 
variance, i.e. the Mott insulating nature of the undoped 
state. Note that the vicinity to the Mott state is also 
responsible for the small supcrffuid stiffness p s ~ x. 

The idea that vortices in the underdoped system have 
SF cores ties in with a more general picture of the pseu- 
dogap state. This is derived from the SU(2) gauge theory 
for the t-J model and views the pseudogap regime as a 
thermally disordered state, where the system fluctuates 
between various types of short range order correspond- 
ing to mean-field states that would all become identical 
at zero doping. The two most prominent correlations are 
d-wave superconductivity - determining the ground state 
as soon as the other fluctuations freeze out at low T - 
and SF correlations. The latter represent - in addition to 
conventional phase fluctuations of the SC order param- 
eter - the lowest lying fluctuations around the d-wave 
SC state with the largest spectral weight 0. We have 
argued^ that the scattering of quasiparticles with these 
SF fluctuations may be related to the partial loss of the 
quasiparticle peaks in the pseudogap state. 

Some of the ideas presented here carry over to other 
types of cheap vortices. Indeed there are experimen- 
tal indicationsTOj for antiferromagnetic ordering at low 
temperatures in the vortex cores of optimally doped Tl- 
compounds, and it is quite likely that the vortex core may 
contain both SF and AF correlations [Tlj. The focus on 
SC and SF correlations is an attempt to concentrate on 
the main tendencies suggested from the SU(2) approach. 
The SF state has the advantage that it has a gap struc- 
ture identical to the d-wave superconductor and naturally 
explains the gap in the vortex core. Other types of cor- 
relations, such as antiferromagnetic tendencies, may be 
viewed as additional instabilities, which naturally coexist 
with the SF order for small doping. 

Recently Ivanov and Lee 01 calculated the energy dif- 
ferences between SF and d-wave SC states using the 
Gutzwiller projection technique. Together with the com- 
puted superfluid stiffness this provides an estimate for 
the energy of a SF vortex. Note that in Ref. ^| the 
vortex core turns out to very small, but on the other 
hand the core state, taken to be a pure SF state, is likely 



to be somewhat too high in energy, and a better core 
state will increase the size of the core again. Keeping in 
mind these uncertainties and for the lack of better pa- 
rameters we will use the numbers of Ref. |l2j as input 
for a generalized XY-model. 

Let us begin with the SU(2) mean-field (MF) theory @. 
Here we are interested in low temperatures. Hence we 
assume that the bosons carrying the electronic charge 
are condensed. The Hamiltonian for the fermionic spin 
degrees of freedom fa, fa on the lattice sites i is given 
by 

X contains the hopping contributions and the con- 
straint field ao,3, which is a function of the doping. 
Next we allow for local admixtures of SF amplitude 
in exchange for SC pairing amplitude, described by 
an angle 0^8], and fluctuations of the SC phase, a«. 
Then the SF amplitude on the bond ij is given by 
Wij — «A (— l)' l * + jy cos 9 ' 1 2 6j and the pairing amplitude 

is A i3 = (-l)vHyA sin^^exp[i(aj +<*j)/ 2 ]- The 
pure superconductor has 9 — n/2, while the two degen- 
erate SF states have 9 — and 9 = it. 

Now consider an effective Hamiltonian for on and 9i, 

H = p s (x, 9i, 9j) cos(a?i — aj) + mg cos 2 9i 

{ij) i^vp 

+Kj2(Oi-0j) 2 + J2 H v- ( 2 ) 

(ij) VP 

The first term is the phase stiffness of the SC phase that 
depends on doping x and the local 9. The second term 
is a mass term for 9 that takes into account the energy 
difference between SF and SC state outside the vortex 
plaquettes (VP). The third term is a gradient term for 
the 0-variation. The last term is the vortex core energy, 
ify = X)igvp( m » sm6 ^ + m e)- The sum is over the four 
sites on each of the VP. Here A^ is assumed to vanish 
and according to Eq. ^ 9 now describes an interpolation 
between the SF state (9 = 0,7r) and the zero flux state 
(9 = it/ 2) which is a Fermi liquid. As discussed earlier, 
the Fermi liquid core is expected to be costly and the 
energy costs of the SF state and the Fermi liquid state are 
mg and mg + m n , respectively. The specific ^-dependence 
chosen for Hy comes from the MF theory for uniform 
9 and 9 is the average 9i over the four sites of the VP. 
Requiring that the vortex core area 4a 2 equals the vortex 
size 7r£ 2 of a pure SF vortex in the microscopic calculation 
in Ref. we obtain that the lattice constant of Eq. 

121 is slightly larger than the one of the underlying t-J 
model with a scale factor of ~ 1.2. We assume that 
the variation of the vortex core size can be neglected 
in the small doping and low temperature region we are 
interested in. 

In principle, all parameters for can be obtained 
from the MF theory by integrating out the fermions as in 
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FIG. 1: Snapshot of the simulation for x = 0.06 and T = 
0.09J « 2T,f. Left: Vortices (circles) and phase angles (ar- 
rows). The arrow length denotes the local SC amplitude. 
The shaded squares indicate sites with large SF admixture, 
| cos Si | > 0.9. The average value of | cos#| in the vortex cores 
is 0.7 (1 for the pure SF state). Right: SF fluctuations cos ft 
of the same sample. 



Rof. 8], but where possible we use parameters obtained 
by the Gutzwiller variational treatment of Ivanov and 
Leepjj. With t/J = 3, this gives p s (x,9 = tt/2) « 
0.75zJ, ra B « 0.33xJ and m„ « (0.25 - x) J/2. The 
^-dependence of p s is obtained numerically from MF 
theory which shows a rapid linear increase from zero 
when 9 deviates from or ir. The coefficient K for the 
gradient of 9 is difficult to extract from the collective 
mode spectrum^, as the #-mode energy is not simply 
quadratic around (%, 7r). Nevertheless the g-dependence 
of the 9- variations is quite weak and less than that of the 
SC phase a and goes to zero for x — > 0. Thus we estimate 
K r~j p s /2. Other choices give similar results. 

In the scaling theory for the BKT transitional the 
vortex core energy is assumed to be large compared to 
the temperature such that the fugacity y = exp{—E c /T) 
can be used as small parameter. In the limit y — * 
the transition occurs at T c = irp s /2. A nonzero fugacity 
y > leads to a reduction of T c from this upper bound by 
roughly ir 2 p s y. With the parameters above the core en- 
ergy of a single ideal SF vortex is E c — Q.75ttxJ or twice 
the maximal T™ ax = np s /2. This leads to a reduction of 
the critical temperature down to T c w 1.06p s . 

The model described by Eq. [21 can be simulated with 
Monte Carlo methods. To estimate T c we calculate the 
helicity modulus T which [l4| measures the rigidity with 
respect to a phase gradient in the system. T vanishes 
above the SC transition: in BKT theory it jumps to zero 
at T c , the height of the jump being 2T c /n, independent 
of the core energy [15j. We also measure the average 9- 
variation inside and outside the vortex cores and how the 
number of vortices depends on doping and temperature. 

A snapshot from the Monte Carlo is shown in Fig. 
for a sample above T c . We observe two relatively well iso- 
lated and other less separated vortices. The SC phase on 
is disordered but exhibits remnants of short range order. 
For T = 0.09 J and x = 0.06, 9 varies rapidly in space due 
to its light mass mg < T and small gradient terms. Hence 
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FIG. 2: Left: Helicity modulus T versus temperature T for 
dopings x = 0.06 and x = 0.12, averaged over 300 100x100 
samples. The temperature where T goes to zero is an upper 
bound for the SC transition temperature T c . The dashed line 
has the slope 2/tt. Right: Vortex density log 10 ny per site vs. 
T and x. The dotted line is T = p s = 0.75zJ. 



there is a significant amount of SF admixture reducing 
the SC pairing amplitude locally even outside the vortex 
cores (see shaded areas in Fig. QJ. Outside the vortices, 
the average SF amplitude is (|cos0|) ~ 0.3 (compared 
to 1 for the pure SF state). Inside the vortex cores it is 
strongly enhanced, but not maximal ((|cos0|) ~ 0.7 for 
the sample shown, see also left plot in Fig. [SJ. 

In the left plot of Fig. [2] we show the helicity modu- 
lus T for two dopings x and a 100x100 system. T goes 
to zero above a doping-dependent temperature, but does 
not exhibit the universal jump of the XY-modelHl E3. 
This is clearly a finite-size effect which is exacerbated by 
the small vortex density and which does not occur in the 
XY-model where the core energy is zero. An estimate for 
the true T c is the intersection of the data with the line 
2T/7T, that is based on the the jump criterion (AY) T = 
7r/2[Tg. With that we arrive at T c e ~ 0.75a;J - p^ A 
numerical bound is T c < xJ . Note that in a <i-wave su- 
perconductor thermally excited nodal quasiparticles lead 
to an additional reduction of the superfluid density that 
is not included here. We expect however that this does 
not affect the nature of the SC transition. 

The temperature and doping dependence of the vortex 
density is summarized in the right plot of Fig. [21 The 
onset temperature for a finite vortex density is an increas- 
ing function of x, approximately given by Ty ~ p s . In 
the underdoped system there is wider temperature range, 
starting at the BKT T c and extending up to the mean- 
field T c ~ 0.18 J where the vortex density continues to 
increase and does not saturate. Thus the vortices do not 
overlap and there is the possibility that SC phase coher- 
ence is still well defined locally in a range above T c . For 
x = 0.06 the phase correlation length £ Q is ~ 7 lattice 
spacings a at T ~ 2T c e . This is in contrast with the nor- 
mal XY-model where at 2T C , £ Q < a, and on average 
there is a vortex on every 5th site. We emphasize that in 
our model the phase fluctuations are not the only source 
of disorder. This can be clearly seen in the right panel of 
Fig. ^ The SF fluctuations are not limited to the vortex 
cores and lead to sizable amplitude fluctuations as well. 
Thus the proximity to the Mott state gives rise to the 
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FIG. 3: Left: SF admixture inside (dashed line) and outside 
(solid line) the vortices vs. temperature T for x — 0.06. Right: 
Density of states for x = 0.06 and T = 0.01, 0.05, and 0.13J, 
averaged over 25 samples (sizes 40x40 up to 48x48). 

low energy scale for SF fluctuations and is responsible 
for both phase and amplitude fluctuations proliferating 
at comparable temperature scales. 

We now address the electronic excitations. At finite 
doping the bulk SF state has small Fermi pockets, and 
correspondingly it is natural to expect that a state which 
fluctuates between a superconductor with gap nodes and 
a SF state with small Fermi arcs will produce a finite den- 
sity of states (DOS) at low energies. When we restrict 
the considerations to static configurations, we can calcu- 
late the quasiparticle spectrum from Hamiltonian for 
a given configuration of fluctuations on a finite system 
and average over many samples. 

Results are shown in the right plot of Fig. [3J The 
DOS exhibits a suppression for all temperatures below 
the mean-field transition at T ~ 0.18J, but the gap fills 
in when T is increased through the SC transition. The lo- 
cal DOS is inhomogeneous, but is not simply correlated 
with or confined to the positions of vortices or regions 
of higher SF amplitude. These results share some as- 
pects with a recent work by Eckl et al. |16( who considered 
disordering the d-wave superconductor by phase fluctu- 
ations modeled by a simple XY model. In their model 
the vortex cores are conventional and the core energy 
is zero by construction. Our generalized model accomo- 
dates both cheap vortices with SF core and energetically 



more expensive vortices with Fermi liquid core. Our cal- 
culation shows that already at the BKT T c there is a fi- 
nite number of quasiparticle excitations at low energies, 
and it is likely that the conductivity will be dominated by 
these normal excitations. The quasiparticle contribution 
to the Nernst signal however is small 1] and the vortex 
contribution will dominate the Nernst signal at low T. 

In conclusion, we have presented a simple model, mo- 
tivated by the SU(2) approach for the high-T c cuprates, 
that describes the superconducting transition as an un- 
binding transition of vortices with staggered flux core. 
Using parameter estimates from projected wave functions 
and the SU(2) mean-field theory it allows us to under- 
stand the occurrence of a moderate vortex density even 
for underdoped systems with low T c . The vortex density 
is determined by an energy scale that is closely related 
to the energy difference between the SF and the d-wave 
SC state and that disappears towards zero doping, mak- 
ing the vortices relatively cheap. However there is wider 
temperature range above T c where the vortices are suf- 
ficiently dilute and do not overlap such that we expect 
that the superconducting phase coherence stays intact 
locally. This is in accordance with the interpretations of 
the Nernst effect measurements 0. Phase and SF (i.e. 
amplitude) fluctuations lead to a filling in of the gap in 
the density of states at small energies already near the 
superconducting T c . This could account for the normal- 
looking in-plane transport of the pseudogap phase which 
onsets just above the superconducting transition. A full 
analysis of this issue requires a calculation of the flux- 
flow resistivity in the fluctuating state. Although our ap- 
proach involves rather strong simplifications and approx- 
imations we believe that it describes a way to understand 
the simultaneous occurrence of normal transport behav- 
ior (e.g. in the resistivity) and strong SC fluctuations, as 
witnessed by the vortex Nernst effect. 
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